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Adiabatic quantum computing and optimization have garnered much attention recently as possible models 
for achieving a quantum advantage over classical approaches to optimization and other special purpose com¬ 
putations. Both techniques are probabilistic in nature and the minimum gap between the ground state and first 
excited state of the system during evolution is a major factor in determining the success probability. In this work 
we investigate a strategy for increasing the minimum gap and success probability by introducing intermediate 
Hamiltonians that modify the evolution path between initial and final Hamiltonians. We focus on an optimiza¬ 
tion problem relevant to recent hardware implementations and present numerical evidence for the existence of a 
purely local intermediate Hamiltonian that achieve the optimum performance in terms of pushing the minimum 
gap to one of the end points of the evolution. As a part of this study we develop a convex optimization formu¬ 
lation of the search for optimal adiabatic schedules that makes this computation more tractable, and which may 
be of independent interest. We further study the effectiveness of random intermediate Hamiltonians on the min¬ 
imum gap and success probability, and empirically find that random Hamiltonians have a significant probability 
of increasing the success probability, but only by a modest amount. 


I. INTRODUCTION 


Controlled adiabatic evolution of a many-body quantum 
system is the basic ingredient for performing quantum com¬ 
puting via adiabatic quantum computing (AQC) 0 or opti¬ 
mization via adiabatic quantum optimization (AQO, also re¬ 
ferred to as quantum annealing if finite temperature effects 
cannot be ignored) I®. The conventional approach to AQC 
and AQO implements a linear interpolation between an initial 
(//q) and final ( I! \) Hamiltonian acting on n spins (or qubits): 

H(t) = (l - Ho + 2jJi, (1) 

where T is the total time of evolution. For convenience, we 
hereby reparametrize the temporal evolution by introducing a 
normalized time parameter s = t/T (0 < s < 1). The system 
begins in the ground state of and the ground state of ll \ 
encodes the solution to a problem of interest. Then evolving 
the system according to H(t) until t = T and measuring per¬ 
mits sampling from the solution state if the system remains 
in the ground state throughout the evolution. The probability 
of this happening, or the success probability of the adiabatic 
computation, relies on satisfying the following adiabatic con¬ 
dition (h= 1 throughout this work) Ob: 

maxo<Ki I (<Ms)| \M S )) I ,,,, 

min 0 < s <iA 2 ( S ) ’ 


Ai(s) > Ao(s), and A(s) = Ai(s)—Ao(s). The quantity A(s) 
is referred to as the instantaneous gap and its minimum value 
over the entire evolution as the minimum gap for the problem. 
Satisfying condition ([2} guarantees, with high probability, that 
the state of the system at the terminal time s = 1 ( i.e ., t = T) 
will be the ground state of the problem Hamiltonian Hi IH- 
It should be noted that the condition in Eq. ({3 comes from a 
worst-case analysis. An alternative, local adiabatic condition 
that arises directly from the adiabatic theorem is HU: 
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( 3 ) 


where A n o(t) = A n (t) — Ao (t). The sum in this condition is 
typically dominated by the n = 1 term. 

While the adiabatic model of quantum computation can 
be shown to be universal lIs UToll . the formulation of a fault- 
tolerant version of AQC remains a challenging and open prob¬ 
lem GMI. Despite this deficiency, there have been sub¬ 
stantial efforts to implement special purpose processors for 
optimization 0,11, and several heuristics for increasing the 
reliability of AQC and AQO have been developed, including 
encoding or increasing the degeneracy of the problem Hamil¬ 
tonian flS - [l8ll . 

A heuristic performance enhancement strategy that has 
been extensively studied is the modification of the interpo¬ 
lation between Hq and Hi from linear as in Eq. 0 to a more 
general time-dependent form: 


where |(^>o(s)), |</>i(s)) are respectively the instantaneous 
ground state and first excited state of H(s), and A(s) is the 
difference in energy between the first excited state and ground 
state at time s, i.e., H(s) |<fo|i( s )) = A 0 |i(s) |0o|i( s )) with 
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H(s) = f 0 (s)H 0 + fi(s)Hi, 

where fo and fi are arbitrary but smoothly varying functions 
satisfying boundary conditions /o(0) = 1, /o(l) = 0 and 
/i(0) = 0, /i(l) = 1. Several authors have shown that 
by tuning the functOions /o|i(s) one can increase the success 
probability of the adiabatic algorithm for a fixed T [fd. IT9l-t24ll . 
The improvement afforded by such a modification of the in¬ 
terpolation scheme can be understood from the intuition that 
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it can be tuned to move slower near smaller energy gaps and 
faster near larger ones, provided that these regions are known 
a priori. Another class of strategies, commonly known as 
shortcuts to adiabaticity j25l 126 1. attempts to optimize fi(s) 
and T to obtain the solution state (ground state of Hi) in the 
shortest time possible, but forfeits adiabaticity during the evo¬ 
lution. 

In this work we study another heuristic strategy for increas¬ 
ing the reliability of AQC and AQO, which involves a further 
generalization of the adiabatic interpolation by adding addi¬ 
tional intermediate Hamiltonians between the two endpoints. 
This leads to a Hamiltonian of the form: 

M+l 

H(s-, {/_,■(*)}) = f 0 (s)H 0 + J2 + h{s)Hi. (4) 

j=2 

The fj are smoothly varying functions of s with /o(s) and 
fi (s) satisfying the same boundary conditions as above, and 
the other M schedule functions satisfying fj(0) = fj( 1) = 0 
for all j > 2. The addition of these extra Hamiltonians pre¬ 
serves the problem structure but introduces the possibility of 
modifying the eigenvalue distributions (and thus gaps) dur¬ 
ing time evolution. In effect, the intermediate Hamiltoni¬ 
ans allow the adiabatic evolution to take paths other than a 
straight line from Hq to Hi. We refer to the collection of 
functions {fj(s)} as the adiabatic schedule and this strategy 
as schedule path optimization (SPO). We will interchangeably 
call H p (s) = intermediate Hamiltonian 

or the (adiabatically) perturbing Hamiltonian since it perturbs 
the straight line path. We note that early attempts to study the 
effects of such intermediate Hamiltonians were made in Ref. 
[27], and more recent investigations are in Refs. lHHE3,[Hf] . 

We will study these intermediate Hamiltonians from the 
perspective of optimization and investigate whether there ex¬ 
ist local intermediate Hamiltonians that increase the success 
probability of computations. We will consider only local, or 
single-qubit, intermediate Hamiltonians for two reasons: (i) 
these are likely the most experimentally feasible; and (ii) the 
number of intermediate terms is restricted to M = 3n at most, 
for an n-qubit problem, and thus the search for optimal per¬ 
turbing Hamiltonians is tractable. We show that the problem 
of finding good intermediate Hamiltonians can be approxi¬ 
mated by a convex optimization, which considerably increases 
the solution efficiency. However, even with a convex opti¬ 
mization approach, explicit solution of the optimal schedule 
path quickly becomes infeasible as n increases. Therefore, we 
also study the likelihood of random local intermediate Hamil¬ 
tonians enhancing the success of computations. In addition 
to AQC and AQO, the approach considered can inform other 
quantum adiabatic protocols such as adiabatic state and pop¬ 
ulation transfer 12911 . In such cases, one may be interested 
in systems with moderate state space dimension, and the ap¬ 
proach we outline can be used to constructively find interme¬ 
diate Hamiltonians that increase the probability of success. 

The remainder of the paper is organized as follows. Section 
Ullintroduces the class of adiabatic evolution problems we will 
focus on, and explicitly formulates the optimization problem 
of interest. Then in section [III] we apply a convex approx¬ 


imation to the optimization of the intermediate Hamiltonian, 
present solutions to the problem, and study their properties. In 
section[IV]we consider the strategy of introducing random in¬ 
termediate Hamiltonians and study its effectiveness. Finally, 
sec tion[Vlsumm arizes the findings and discusses implications. 

II. PRELIMINARIES 

For our investigations we focus on adiabatic interpola¬ 
tions that solve quadratic unconstrained binary optimization 
(QUBO) problems. This NP-hard problem [[30] is phrased as 
a minimization of a sum clause of binary variables, for exam¬ 
ple: 

min axi + bx 2 + CX 1 X 2 
£ 1,212 

s. t. xi, x 2 £ {0,1}. 

A quantum approach to solve such a problem with n vari¬ 
ables is to map the solution to the ground state of the final 
Hamiltonian of an n-qubit system, and then let the quantum 
system adiabatically evolve from the ground state of a readily 
prepared initial Hamiltonian to the ground state of that final 
Hamiltonian. In particular, we can choose the final Hamilto¬ 
nian as 

Hi =Y^ h i a l + *3*1*1, (5) 

i i,j 

where the eigenstates of a\ correspond to the binary values of 
variable Xi, and the Hamiltonian coefficients hj and J, :l are 
functions of the problem coefficients (a, 6, c in the example 
above). Note that we suppress tensor product symbols and 
identities when denoting n-qubit Pauli matrices, and super¬ 
scripts indicate the qubits on which there is a non-trivial ac¬ 
tion; i.e., a%. = I <g> a x <8> I ® • ■ ■ <g> I. The initial Hamiltonian 
can be set in the standard form 

= ( 6 ) 

i 

Although it is not possible to transform any classical computa¬ 
tion to a QUBO problem these problems are particularly rele¬ 
vant since recently developed experimental devices attempt to 
solve problems from this class [0, |5|] . 

A. Problem specification 

In order to assess the effectiveness of local intermediate 
Hamiltonians, we must define a performance metric. One 
would ideally use the probability of success p sncc (T ) = 
\{(f>o(T)\tjj(T))\ 2 , where |</>o(T')) is the ground state of Hi 
and \tjj(T)) is the actual achieved state at the final time. How¬ 
ever, for the ease of optimization formulation, we instead use 
a proxy for the success probability, namely, the minimum gap 
A m in = min , A(s). This enables us to approximate the inter¬ 
mediate Hamiltonian optimization problem by a convex op¬ 
timization, which facilitates numerical investigation of inter¬ 
mediate size qubit systems. In Sec. QTi]we will confirm that 
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the optimal intermediate Hamiltonians identified by this proxy 
cost function do indeed enhance the probability of success for 
large T in most cases. 

Now we are in a position to formulate the optimization 
problem. Given an adiabatic problem specification with 
Hq, Hi, and a set of M local, intermediate Hamiltonians 
{Hj}^^ 1 , the optimal schedule path is defined as the one 
that solves: 

u^m 1 °<Si A(s;{/j(s)}) 

s.t. /j(0) = /j(l) = 0, Vj> 2, (7) 

\fj(s)\<f B , Vj >2, 

l/i(s)|<£, V j > 2, 

where A(s; {/j(s)}) is the gap of fT(s; {fj(s)}). In other 
words, the optimal schedule path maximizes, over the set of 
interpolation functions, the minimum gap on the entire time 
duration. The first set of constraints forces the intermediate 
Hamiltonians to be zero at the initial and terminal time, the 
second set bounds the amplitude of the intermediate terms, 
and the last set ensures that the control functions are smooth (e 
is a small positive constant). We fix/o and/i as/o(s) = 1 — s 
and /i(s) = s for simplicity and also to isolate the effects 
of the intermediate Hamiltonians from effects of varying the 
interpolation velocity. 

The minimax optimization above is difficult to tackle in its 
continuous form. Therefore we first discretize time to obtain 
a more tractable formulation. Divide the total time duration 
s £ [0,1] into N equal intervals each of length As = A 
Within the interval [*As, (i + l)As), assume that the sched¬ 
ule functions fj(s) take a constant value fj(i) = fj(iAs), 
where * = 0, • ■ ■, N — 1. In particular, / 7 (0) = fj(N ) = 0. 
When N is large enough, this piecewise constant function ap¬ 
proximates the continuous function well. We then define a 
M(N — 1) dimensional discrete optimization variable 

A = [/ 2 (1), /a(2), f m+i{N - 2), f M+ i(N - 1)] 

and the discrete optimization problem can be written as 

max min AAA) 

A l<i<N-l 

S. t. \fj{i)\ < f B , i = !,••• ,N- 1, 

I fj(i + 1) - < eA s, i = 0, • • • ,N - 1, 

where j = 2, • • •, M + 1. The cost function A i(A) is the 
energy gap between the ground and the first excited states at 
s = iAs when the particular adiabatic schedule defined by A 
is used. 

For small number of qubits, this discrete optimization can 
be solved efficiently using standard sequential quadratic pro¬ 
gramming methods. In particular, we applied the Matlab Op¬ 
timization Toolbox to solve for the (not necessarily globally) 
optimal intermediate Hamiltonian for a 10-qubit QUBO prob¬ 
lem in 82.9 hours on an ordinary desktop computer. However, 
when the number of qubits increases further, the numerical 
solution quickly exceeds the computational power of avail¬ 
able computers. To remedy it, we will formulate a convex 
optimization approximation to Eq. © in the next section. 


Before proceeding we note that this method of gap opti¬ 
mization cannot be utilized to increase success probability if 
the minimum gap for the problem occurs at the endpoints of 
the interpolation (i.e., at s = 0 or 1). However, most QUBO 
problems do not have minimum gaps at s = 0 or 1 under lin¬ 
ear interpolation. 


III. OPTIMAL INTERMEDIATE HAMILTONIANS FROM 
CONVEX OPTIMIZATION 


We now approximate the search for optimal adiabatic 
schedules by a convex optimization formulation motivated by 
the methods in Ref. 113ltl . Convexity ensures that any local 
optimum is a global optimum and a convex problem can be 
solved efficiently with mature numerical methods j!32J|. Thus 
this approximation technique makes it possible to numerically 
search for optimal adiabatic schedules for larger AQC prob¬ 
lems; e.g., QUBO problems up to 22 qubits readily. 


A. Approximate convex optimization formulation 


Denote the Hamiltonian at time instant i As as It (i), and its 
eigenvalues and eigenvectors as Afc(z) and u k (i), i-e., 

H(i)u k {i) = A k {i)u k {i), i = 0, • • • , N. 


For brevity of notation, we drop the dependency on the op¬ 
timization variable A. Since H(i) is Hermitian, we assume 
that Uk(i) are chosen as an orthonormal basis and A k(i) are 
arranged in ascending order. Consider an upper bound eo(i) 
for Ao(z) and a lower bound ei(z) for Ai(z). The discrete op¬ 
timization problem Eq. © can be written as 


,iV -1, 


s. t. A 0 (z) < e 0 (z), i = l,-" 

Ai(z) > ei(z), i = 1,••• ,N — 1, 
\f,(i)\<f B , i = N-l, 


(9) 


I h 


\ / | J 1 - — 7 7 — 7 

(* + 1) - fj (*)I < eAs, i = 0, • • - ,N- 1, 


where j = 2, • • •, M + 1. The first two constraints in ©. 
i.e., the bounds on Aq (i) and Ai(z), can also be written as 


$o(*) (#(*) - e o(*)I) ®o(i) < 0, 

$i(*) (H(i) - ei(z)I) ^i(i) >r 0, 

where 


( 10 ) 


$o(*) = [uo(*)], (11) 

$i(z) = [tii(z), ■ •• ,U2»-i(z)]. (12) 

Here A is the Lowner partial ordering on positive semidefi- 
nite matrices, i.e., A A B if and only if A — B is positive 
semidefinite. 

We can use Eq. (fTTTb to substitute the first two inequality 
constraints in © and obtain an equivalent optimization prob¬ 
lem. This is still a nonconvex problem, because the con¬ 
straints ([Tol l are nonconvex. However, we can modify the 
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formulation slightly by using an iterative procedure where the 
eigenvectors Uk{i) from the previous iteration are used to form 
the matrices $o(*) and $i(i) for the current iteration. This 
way, $o(*) and $i(*) become constant matrices and do not 
depend on the current Hamiltonian. We then obtain a convex 
optimization problem since both the cost function and con¬ 
strains are convex. Provided that the iteration step size is small 
enough, it can well approximate the original adiabatic sched¬ 
ule search problem (more about this approximation later in 
this section). 

When dealing with a system with n qubits, the matrix $1 (?) 
has 2" — 1 columns, which may demand excessive computa¬ 
tional power for large value of n. We thus drop the eigenvec¬ 
tors in T] (?') that correspond to higher energies to reduce the 
problem size m, that is, we truncate 

,u p (i)], p<2" —1. (13) 

In our numerical studies of systems up to n = 22, we have 
observed that p = 5 suffices to yield good results. 

Now after these two approximations, we have the following 
reduced approximative convex optimization problem formula¬ 
tion: 


max min ei (i) 

A,e 0 (i),ei(j) l<i<JV-l 


e 0 (*) 


s. t. $J(i) (H(i) - eo(i)I) fcoW < 0, * = 1, • ■ ■ , N - 1, 
$!(*) (H(i) - ei(i)I) $i(«) ^0, i = 1, • • • , N - 1, 
< / b , * = ,N- 1 , 

\fj(i + !) ~ fj(i )I < eA s, i = 0, ,7V- 1, 

(14) 

where is dehned in Eq. (fill and <!>{(*) in Eq. (fl3l) . The 
complete algorithm to solve for the optimal adiabatic schedule 
can be summarized as follows: 


Step 1 Set the initial guess A = 0 and let A = A in the first 
iteration. Set the values of two small positive constants 
£ and r). Calculate the minimum gap Ao with linear 
interpolation. 

Step 2 Compute the matrices and $i(i;-A). Solve 

the convex optimization problem (fl4l > with the addi¬ 
tional convex constraint ||A — A|| < rj to ensure the 
resulting optimal solution A* close to A. 


Step 3 If II min A (i: A*) — min A(i; A)|| < £, stop the 

l<i<N -1 i=0 ,n ~ 

algorithm and return the optimal solution A*; other¬ 
wise, let A = A* and go back to Step 2. 


The additional convex constraint specified in Step 2 of this al¬ 
gorithm aids in restricting the change in the optimization vari¬ 
able A between iterations, which is necessary for the eigen¬ 
vectors Uk(i) of previous iteration to be a good approximation 
to the eigenvectors of the Hamiltonian for the current itera¬ 
tion. This is obviously not a sufficient condition to ensure the 
validity of this approximation all the time, but in practice we 
have found that the above algorithm performs well when this 
condition is imposed. 


B. Optimal adiabatic schedules 

We now present the results of our search for optimal local 
intermediate Hamiltonians for the QUBO problem. We in¬ 
vestigate a variety of QUBO problems, ranging from 3 to 22 
qubits in size. For n < 10 qubits we directly solve the min¬ 
imax optimization problem formulated in Eq. ©; however, 
for larger problem sizes we solve the convex approximation 
developed above. For each problem size, we do not assume 
any special structure for the local terms or the couplings, but 
rather, we generate random problem Hamiltonians by sam¬ 
pling each hi and ./ y in Eq. © from a uniform distribution on 
the interval [—1,1]. In the following, we present results from 
a 19-qubit and a 22-qubit system. These results are represen¬ 
tative of the general properties of solutions to all the problems 
we investigated. 

The initial Hamiltonian is set to be Hq = ]>A a' x , and thus 
the gap at the initial time is 2, which is the cost of flipping 
the phase of any one qubit. For the two examples considered 
here, we choose the final Hamiltonian Hi to also have a gap 
around 2 and the minimum gap under linear interpolation to be 
located away from the two endpoints. These choices are made 
to demonstrate the effect of optimal adiabatic schedules more 
effectively. The intermediate Hamiltonian takes the form 

n 

ifQ UBO ( s ) = ]T {f 2 0 (s)ai + f 2j +i{s)ai) , (15) 

3 = 1 

and hence M = 2 n. We fix the transverse field in the x di¬ 
rection and ignore local a y perturbations in Hq and II p . This 
is because the sum of a local a y and a x term at any s is still 
a net local field in the transverse direction orthogonal to the 
problem Hamiltonian and thus a local tr y term does not have 
any additional effects. 

We also divide the time interval [0,1] into N = 50 equal 
steps, and set the maximum allowed changing rate of sched¬ 
ule functions as e = 2.5. Numerical studies reveal that vary¬ 
ing the discretization number N of the time interval has little 
effect on optimization results as long as its reasonably large, 
and we shall explore the effect of varying e later. Finally, the 
maximum value of the intermediate terms fn is set to be the 
maximal value of \hi \ and | |, ensuring that the energy of the 

local terms in the intermediate Hamiltonian is never excessive 
compared to the terms in the problem Hamiltonian. 

The results of solving Eq. i fTTl i for a 19-qubit QUBO in¬ 
stance are shown in Fig. Q] After 2 iterations and 2.2 hours 
on a single processor we obtain adiabatic schedules that can 
push the minimum gap to the terminal time. Fig.[2]shows the 
results for a 22-qubit instance after 5 iterations and 112 hours. 
In this case, the minimum gap increases with each iteration 
of the convex optimization as the adiabatic schedule is refined 
further. Most of the computing time in both cases is devoted 
to calculating eigenvectors and eigenvalues as opposed to per¬ 
forming the convex optimization. 

A striking feature present in both solutions in Figs. Q]and 
[2]is that for the (locally) optimal adiabatic schedule, the min¬ 
imum gap is equal to, or almost equal to, the minimum of 
the initial and final gaps. We refer to this as the best-case 
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FIG. 1: (Color online) Optimized adiabatic schedules for a 19-qubit 
AQC problem. Top: Energy gap, A(s) for linear interpolation ac¬ 
cording to Eq. CD (blue dashed), and after SPO (black solid). Bot¬ 
tom: The M = 38 schedule functions, for the optimized adi¬ 

abatic schedule. The dashed lines show the time-dependent coeffi¬ 
cients of local a x terms and the solid lines show the coefficients of 
local a z terms. 



Time (s) 


FIG. 2: (Color online) Optimized adiabatic schedules for a 22-qubit 
AQC problem. Top: Energy gap, A(s) for linear interpolation ac¬ 
cording to Eq. 0} (blue dashed), and after SPO (black solid). Bot¬ 
tom: The M = 38 schedule functions, fj(s ), for the optimized adi¬ 
abatic schedule. The dashed lines show the time-dependent coeffi¬ 
cients of local a x terms and the solid lines show the coefficients of 
local <r z terms. 


performance (i.e., min s A(s) = min{A(0), A(l)}) since the 
minimum gap has been pushed to one of the endpoints, where 
H p = 0. This is seen to be a generic feature of the optimal 
adiabatic schedules for all the QUBO problems we studied. 
In the non-convex formulation, which is solvable for n ~ 10 
qubits, the optimal adiabatic schedule achieved this best-case 
performance for all the instances that we solved. For the con¬ 
vex approximation, which is solvable for larger n, this best- 
case performance is approached with increasing iterations. 

We have not been able to discern any pattern in the form of 
the optimal schedule functions other than that some of them 
generally peak near the region of minimum gap for the linear 
interpolation and decay (mostly smoothly) to zero, whereas 
some others go across zero near this region. 


C. Effect of optimization parameters 

The primary free parameter in all formulations of the gap 
optimization (Eqs. ([?}, ©, and ( 1 1 -H i) is the bound on the rate 
of change of the schedule functions, e. Allowing for large 
slew rates increases dH / ds, which in turn increases the run¬ 
time of the adiabatic algorithm, as illustrated by Eq. ©. At 
the same time, the larger this parameter is the more freedom 
there is in varying adiabatic schedule functions, and this can 


1.6 
1.4 
1.2 
1 

0.8 

0 0.5 1 0 0.5 1 1.5 2 2.5 

Time (s) e 

FIG. 3: The effect of e on the optimization performance for the 19- 
qubit example discussed in Fig. Q] Left: the gap as a function of 
s for linear interpolation; Right: the minimum gap achieved by the 
optimization as a function of e. 


enable greater tuning of the system minimum gap. 

This is evidenced in Fig. [3] where we show the mini¬ 
mum gap achieved by SPO for the 19-qubit problem pre¬ 
sented in Fig. Q] as a function of e. The optimization is per¬ 
formed using the convex optimization formulation ( fl4l >. We 
see that for small e the maximum minimum gap achieved by 
the optimization is modest, and then as the value of £ in¬ 
creases the minimum gap saturates to the best-case perfor¬ 
mance {i.e., min s A(s) = min{A(0), A(l)}). However, the 
minimum value of e to achieve this best-case performance is 
problem dependent and cannot be determined a priori. 
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D. Relation to AQC success probability 


We have demonstrated that local intermediate Hamiltoni¬ 
ans can increase the minimum gap of QUBO problems; how¬ 
ever, as mentioned in the Introduction, the minimum gap is a 
proxy for the success probability of an AQC. Increasing the 
minimum gap is likely to increase success probability (given 
a fixed running time T), but the relationship between these 
quantities has been shown to be non-monotonic in some cases. 


e-g.. 


We now examine the relation between the success probabil¬ 
ity for the optimized and linear interpolation adiabatic sched¬ 
ules. Tab le|T] I ists the success probabilities for several problem 
instances and total evolution times, with the optimized adi¬ 
abatic schedule as determined by maximizing the minimum 
gap. We see that in most cases, increasing the minimum gap 
also increases the probability of success, with one exception 
when the running time is short (case 1, T = 5). We deter¬ 
mined that this exception is an example of where the increase 
in minimum gap is not sufficient to compensate for the in¬ 
creased variation produced by SPO. That is, the optimized adi¬ 
abatic schedules typically have greater variation in the Hamil¬ 
tonian, | dH/ds |, than the linear interpolation, and for short T, 
the increase in the minimum gap may not be enough to com¬ 
pensate for this increase in variation. This is the reason that 
as T is increased the success probability under SPO exceeds 
that of linear interpolation. However, in the other extreme, 
where T tends to infinity, there is also no advantage to SPO 
since for arbitrarily long running times the minimum gap is 
inconsequential (we do not show data for such long running 
times in Table U but have noticed this effect in the numerical 
studies). Of course this regime is impractical since other fac¬ 
tors and timescales usually set a maximum affordable running 
time. 


E. Hard QUBO instances 

In this subsection we study very hard instances of QUBO, 
and whether SPO can improve performance for these cases. 
To do so, we generated 180,000 random instances of 10-qubit 
QUBO problems and selected the 30 instances with the small¬ 
est success probability for linear interpolation between Hq 
and H\ with T = 10. For conciseness, in the following we 
only present data for the three hardest from these 30 instances 
(ones with the lowest success probability under linear inter¬ 
polation), but the conclusions we draw are based on studying 
the large set of all 30 hard instances. The SPO optimizations 
were performed using the convex approximation in Eq. d 
with e = 2.5. 

Fig. Q] shows the time-dependent gap for the three hard¬ 
est instances with linear and SPO interpolations. These three 
cases exhibit two features that are common to all of the 30 
hard instances we investigated. The first feature is that with 
linear interpolation the minimum gap always occurs near s = 
1 for the hard problems. The second feature is that with the 
optimized schedule, the minimum gap can always be pushed 
out to s = 1; or in other words, SPO (with single qubit inter- 



la) Hard case 1: A^ in = 0.0222, s m in ~ 0.86, Pa UCC = 0.0035. 



(b) Hard case 2: A° lin = 0.0006, s m in ~ 0.96, p° ucc = 0.0066. 



(c) Hard case 3: A^ in = 0.0044, s m in ~ 0.96, Ps UCC = 0.0008. 

FIG. 4: (Color online) The time-dependent gap for the three hardest 
instances of QUBO we found, with linear and SPO interpolations. In 
the above, A^, in is the minimum gap under linear interpolation that 
occurs at s m in, and Pg UCC is the probability of success under linear 
interpolation with T = 10. The inset in (b) shows the overlap of 
neighboring ground states, |(uo(*)|wo(* + 1))|. The sudden drop in 
this quantity around s m i n indicates a finite-size analog of a phase 
transition. We do not show this overlap/fidelity for the other problem 
instances since they do not possess a phase transition, and hence do 
not show a significant change in this quantity when crossing s m in- 


mediate Hamiltonians) achieves best-case performance even 
for the hard QUBO instances. 

We note that hard case 2 seems to undergo a finite-size ana¬ 
log of a phase transition around the minimum gap region un¬ 
der linear interpolation. This is indicated by an abrupt change 
in the ground state as s m ; n is crossed, as indicated by the rapid 
decay of fidelity between neighboring (in s ) ground states 
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Success probability 

10-qubit 

12-qubit 

19-qubit 

Case 1 

Case 2 

Case 3 

Case 4 

Case 5 

Case 6 

Case 7 

lineal' 

SPO 

lineal* 

SPO 

linear 

SPO 

lineal* 

SPO 

lineal* 

SPO 

lineal* 

SPO 

linear 

SPO 

T = 5 

0.2441 

0.1838 

0.1660 

0.5418 

0.2846 

0.6856 

0.1341 

0.3182 

0.1398 

0.1861 

0.1846 

0.8748 

0.1011 

0.2567 

T = 10 

0.3608 

0.8746 

0.3338 

0.9223 

0.5857 

0.9761 

0.2738 

0.9553 

0.2013 

0.9203 

0.3584 

0.9604 

0.3898 

0.8097 

T = 20 

0.4259 

0.9885 

0.5810 

0.9843 

0.8194 

0.9903 

0.3907 

0.9942 

0.2583 

0.9702 

0.5822 

0.9872 

0.7506 

0.9888 

T = 40 

0.4714 

0.9947 

0.8282 

0.9981 

0.9633 

0.9959 

0.4561 

0.9984 

0.5054 

0.9964 

0.8152 

0.9950 

0.9619 

0.9957 

Minimum gap 

0.1106 

1.6704 

0.2989 

1.6549 

0.4728 

1.6659 

0.1381 

1.5491 

0.2245 

1.5862 

0.4231 

2.000 

0.4823 

1.6135 

Time min. gap occurs 

s = 0.50 

s = 1 

s = 0.44 

s = 1 

s = 0.38 

s = 1 

s = 0.44 

s = 1 

s = 0.50 

s = 1 

s = 0.38 

s = 0 

s = 0.34 

s = 1 


TABLE I: Success probability for different running times T, with linear and SPO interpolations, for seven problem Hamiltonians. For each 
case, the minimum gaps under linear interpolation and the optimized schedule are shown, as well as the time when the minimum gap occurs. 


133, 351], which is shown in the inset of0Jb). Such problem in¬ 
stances are known to be extremely difficult for AQO l l36l - (39ll . 

Furthermore, in Fig. [5] we show the probability of success 
as a function of T for the three hardest instances when the lin¬ 
ear interpolation or the optimized path is used. As in the easy 
QUBO cases analyzed above, for short T the small increment 
in minimum gap does not directly translate to an increased 
probability of success. However, for large T the optimized 
schedule does have a greater probability of success than linear 
interpolation. Notably, for hard case 2 the success probabil¬ 
ity cannot be significantly increased even for large T\ i. e ., the 
difference in success probability between linear interpolation 
and SPO is not as significant for this case. We can under¬ 
stand this by noting that the phase transition, and associated 
avoided crossing, in this problem instance is not removed by 
the optimized adiabatic schedule (as evidenced in Fig. |4}b) 
by the persistence of fidelity decay even with SPO), and thus 
the problem remains difficult for AQO. Also, we note that for 
hard cases 1 and 3, the success probability under linear inter¬ 
polation can be dramatically increased by increasing the run¬ 
ning time T. However, for fixed running time, especially for 
T ~ 10, these are difficult for AQO. 

We submit the results in this subsection and in sub¬ 
section ME as numerical evidence that min s A(s) = 
min{A(0), A(l)} can always be achieved for QUBO with 
purely local intermediate Hamiltonians of the form Eq. ( IT5 ] i. 



FIG. 5: (Color online) Success probability as a function of running 
time T for the three hardest instances of QUBO we found, with lin¬ 
ear and SPO interpolations. For long T, SPO outperforms linear 
interpolation, even approaching p SU cc = 1 for hard case 1 and 3. 
Note that the a>axis is not linear. The non-monotonic behavior of 
success probability with T exhibited by the above curves is a well 
established phenomenon; non-adiabatic, fast sweeps may have better 
performance for hard problem instances f24l [40[|. 


IV. RANDOMIZED INTERMEDIATE HAMILTONIANS 


The results in the previous section provide strong evi¬ 
dence that there exists at least one local intermediate Hamilto¬ 
nian for every QUBO problem that achieves best-case perfor¬ 
mance, i.e., min s A(s) = min{A(0), A(l)}. However, find¬ 
ing this intermediate Hamiltonian is a difficult task in general. 
Therefore in this section we ask how well randomly chosen in¬ 
termediate Hamiltonians can increase the minimum gap of the 
linear interpolation adiabatic evolution. Crosson et al. consid¬ 
ered a similar question in Ref. (24 ll . where they chose inter¬ 
mediate Hamiltonians of the form 


7Tp Uad (s) 


= s(i-«)Z! 

j =1 Ot=X,Z 



73 
' a. 5 


(16) 


where the coefficients c‘ t are chosen randomly according to 
a Gaussian distribution. We have evaluated several forms 


for the random intermediate Hamiltonian, but will present re¬ 
sults from sampling over intermediate Hamiltonians of the 
same form as Eq. (IT6l) except for two differences; (i) whereas 
Crosson et al. considered a° a that are one- and two-qubit oper¬ 
ators, we will use only local (one-qubit) perturbations, and (ii) 
we will sample c l a over the uniform distribution over [—1,1], 
(0,1], or [—1,0). This form for the intermediate Hamiltonian 
has the benefit that it is smoothly varying as a function of s for 
any random instance, and also the number of random variables 
needed is 2 n, as opposed to scaling with the discretization of 
the s variable. This makes sampling over such random inter¬ 
mediate Hamiltonians much more feasible. 

In this section we restrict our study to 10-qubit problems, 
and evaluate the change in the minimum gap achieved by the 
intermediate Hamiltonian: 

^ — ^miii — A° lin , (17) 

where A|] ljn is the minimum gap for linear interpolation, and 
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A^n is the minimum gap for the Hamiltonian H(s) = (1 — 
s)H 0 + sHi + .Hp* uad (s). We will also investigate the change 
in success probability: A p = pf ucc -p° ucc , where pf ucc (p° ucc ) 
is the success probability using the perturbing Hamiltonian 
(linear interpolation) with T = 10. 


A. Average case behavior 

First, we evaluate the effect of random intermediate Hamil¬ 
tonians on randomly sampled QUBO instances. We gener¬ 
ate N s = 100 instances of QUBO problem Hamiltonians by 
sampling hi and ,/ y from a uniform distribution over [—1,1], 
and for each one of these instances we assess the increase in 
minimum gap achieved by lV r = 5000 random intermediate 
Hamiltonians of the form Eq. (fl~6] >. with c l a sampled over the 
uniform distribution over [—1,1]. Fig. [6] shows the distri¬ 
bution of 0 and Ap succ over this random sampling of prob¬ 
lems instances and perturbing Hamiltonians. If we average 
over the problems instances and the perturbing Hamiltonians, 
the average-case performance is summarized in Table |lT] We 
see that the random perturbing intermediate Hamiltonian of 
the form in Eq. (IT6]> does not increase performance with high 
probability, and on average seems to decrease the performance 
of AQC. 

Next, we modify the character of the random perturbing in¬ 
termediate Hamiltonian by restricting the domain over which 
the coefficients c l a in Eq. (fl~6]> are sampled. The four restric¬ 
tions we consider are: 

1. All d a positive; 

2. All c l a negative; 

3. c l a corresponding to local cr- terms are positive whereas 
the ones corresponding to local a x terms are negative; 

4. c l a corresponding to local er z terms are negative whereas 
the ones corresponding to local a x terms are positive. 

For brevity, we do not present the full distribution of outcomes 
for each of these cases, but rather just summarize the average 
behavior in table[lll](the averages are taken over the same 100 
QUBO problem instances, and 1000 random perturbations in 
each case). We see that the average behavior with restricted 
c l a can be much better than when the domain of c' n is over 
[—1,1]. In particular, when all c l a are sampled from the pos¬ 
itive interval (0,1] the relative frequency of instances where 
the random intermediate Hamiltonian increased the success 
probability is greater than 0.5, although the actual increase in 
success probability is modest. 


B. Hard QUBO instances 



Instance number 
(a) Distribution of f2. 



Instance number 
(b) Distribution of Ap S ucc- 


FIG. 6: (Color online) The change in minimum gap (a) and success 
probability (b) for 100 QUBO problem instances under 5000 random 
perturbations of the form in Eq. COD, when run for T = 10. 


in gap and in success probability, when 1000 random interme¬ 
diate Hamiltonians of the form Eq. dT6l) are used to perturb the 
linear interpolation schedule. We see that the gap increases on 
average and moreover, increases in most cases for all three 
hard cases. Furthermore, the success probability increases 
with very high probability over the 1000 perturbing Hamil¬ 
tonian samples for all three hard cases. This is in agreement 
with the observations in Ref. 0, although their interme¬ 
diate Hamiltonians take a slightly different form. Therefore 
we have evidence that for these hardest instances of QUBO 
the strategy of introducing random local perturbations to the 
schedule path increases the performance with high probabil¬ 
ity, even though the degree of increase (in success probability) 
is very modest, at least for the form of intermediate Hamilto¬ 
nians considered here. 


Now we return to the hardest instances of QUBO that we 
identified in section UHEl and assess the effect of random inter¬ 
mediate Hamiltonians of the form Eq. (fl6t (with c l a sampled 
uniformly from the interval (0,1]) on solving them using adia¬ 
batic evolution. Figs. [7]and[8]show histograms of the changes 


V. DISCUSSION AND CONCLUSIONS 

We studied the feasibility of increasing the success proba¬ 
bility of AQO implementations using intermediate Hamiltoni¬ 
ans to modify the adiabatic schedule. We began by examining 
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Change in minimum Gap, H 

Change in success probability, Ap SUC c 

Mean value and average percentiles 

-0.014828 (-0.0453,0.0195) 

-0.016424 (-0.0542,0.0163) 

#increase 

7 

12 


TABLE II: Change in minimum gap and success probability over 100 QUBO problem instances each with 5000 perturbing Hamiltonians. The 
first row shows the average change (averaged over problem instances and perturbing Hamiltonians) and average 35th and 65th percentiles 
in parentheses. For each problem instance these percentiles are computed over the 5000 perturbing Hamiltonians and then these percentiles 
are averaged over the 100 problem instances (this provides a measure of the variation in performance of random perturbing Hamiltonians, 
averaged over problem instances). The second row shows the number of problem instances (out of 100) where the average gap or average 
success probability increases. 


Ca sampling restriction 

all + 

all 

O x +, <J z '■ — 

Oz • H - , Ox • — 

n 

-2.2 x 10~ 5 

-0.0213 

0.0022 

-0.0214 

(-0.0144,0.0160) 

(-0.0431,-0.0118) 

(-0.0104,0.0191) 

(-0.0435, -0.0124) 

Ap s ucc 

0.017451 

-0.040459 

0.006233 

-0.028102 

(0.0064, 0.0404) 

(-0.0714, -0.0373) 

(-0.0094, 0.0264) 

(-0.0538, -0.0214) 

# of increased minimum gap 

56 

27 

44 

41 

# of increased success probability 

68 

22 

55 

31 


TABLE III: Change in minimum gap (H) and success probability (Ap succ ) over 100 QUBO problem instances each with 1000 perturbing 
Hamiltonians, for each type of perturbing Hamiltonian presented in the main text. The first two rows show the average change (averaged 
over problem instances and perturbing Hamiltonians) and average 35th and 65th percentiles in parentheses. For each problem instance these 
percentiles are computed over the 1000 perturbing Hamiltonians and then these percentiles are averaged over the 100 problem instances (this 
provides a measure of the variation in performance of random perturbing Hamiltonians, averaged over problem instances). The last two rows 
show the number of problem instances (out of 100) where the average gap or average success probability (averaged over the 1000 random 
intermediate Hamiltonians) increases. 



1000 


800 

C 600 
(U 

3 

~ 

0) 400 

£ 

200 

° -- 1 • ■ 0 
Change in minimum gap 

(c) Hard case 3: mean 
7.85 x 10 -4 . 



200 


150 



0 0.05 0.1 

Change in success probability 


(c) Hard case 3: mean 0.0285. 


FIG. 7: (Color online) Histograms of the change in minimum gap 
for the 3 hardest QUBO problems induced by random intermediate 
Hamiltonians. The caption for each subfigure shows the average gap 
change. 


FIG. 8: (Color online) Histograms of the change in success probabil¬ 
ity for the 3 hardest QUBO problems induced by random intermedi¬ 
ate Hamiltonians. The caption for each subfigure shows the average 
change in success probability. 


the minimum gap during the interpolation between Hq and 
Hi as a proxy for success probability, and developed an ap¬ 
proximative convex optimization formulation of the problem 
of finding intermediate Hamiltonians that achieve the maxi¬ 


mum minimum gap, a problem we refer to to schedule path 
optimization (SPO). Then by explicitly optimizing over inter¬ 
mediate Hamiltonians we collected convincing numerical ev¬ 
idence that there exists at least one purely local intermediate 
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Hamiltonian that achieves the best-case performance (where 
the minimum gap is equal to the gap at one of the endpoints) 
for QUBO problem Hamiltonians. Next we evaluated the ef¬ 
fect of SPO on the hardest QUBO problems, which are the 
ones with minimum gap very close to the endpoint (s = 1), 
and concluded that although SPO can achieve an increase in 
gap and success probability, the improvements are modest. 

Although we have evidence that SPO can achieve best-case 
performance, it is obviously not feasible for real AQC or AQO 
implementations. Therefore in section[IV]we focused on ran¬ 
dom local intermediate Hamiltonians, sampled from a partic¬ 
ular quadratic form of the schedule. We showed that if the 
random intermediate Hamiltonian is sampled correctly (with 
positive coefficients for each local term), the performance of 
AQO on randomly sampled QUBO problems can be improved 
with probability greater than 0.5. This strategy also provides a 
very modest improvement in success probability for the hard¬ 
est QUBO problem instances. More work is required to de¬ 


velop a heuristic strategy that improves performance for these 
very hard QUBO instances. 
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